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Abstract 



' It is shown that the equation of state of nuclear matter can be determined 

o ■ 



within the mean-field theory of guj model provided only that the nucleon ef- 
fective mass curve is given. We use a family of the possible nucleon effective 
mass curves that reproduce the empirical saturation point in the calculation of 
the nuclear binding energy curves in order to obtain density-dependent effec- 
tive coupling constants. The resulting density-dependent coupling constants 

x 

5_i ■ may be used to study a possible equation of state of nuclear system at high 

density or neutron matter. Within the constraints used in this paper to M* of 
nuclear matter at saturation point and zero density, neutron matter of large 
incompressibility is strongly bound at high density while soft neutron matter 
is weakly bound at low density. The study also exhibits the importance of 
surface vibration modes in the study of nuclear equation of state. 
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I. INTRODUCTION 



The relativistic aoo model in the mean-field theory has been very successful in describing 
the properties of nuclear matter and finite nuclei [|l]-[J. It provides a useful framework to 
evaluate the one-baryon-loop vacuum fluctuation effects at finite density. In the mean-field 
approximation one can obtain the baryon self-energies in nuclear matter taking into account 
the Dirac sea modified due to the presence of valence nucleons. The energy density contains 
the modified negative energies, the correction that removes half of the interaction energy 
between the states in the Dirac sea, and the interaction energy between the valence particles 
and those in negative energy states. 

Recently, Brockmann and Toki |J reported a relativistic density-dependent Hartree cal- 
culation for finite nuclei, where the results agree very well with experiment. They employ 
a and u exchange for the effective interaction without taking into account the vacuum fluc- 
tuation effects and adjust the density-dependent effective coupling constants such that the 
relativistic Hartree calculation reproduces the nuclear matter Dirac-Brueckner-Hartree-Fock 
(DBHF) results based on a realistic NN interaction. The density- dependent coupling con- 
stants may effectively contain exchange term and contributions of iterated meson exchange. 
The effective scalar field for instance, contains large contribution from two pion exchange 
and the iterated exchange of tt- and p-mesons contribute to the effective vector fields. 

We show that the energy per nucleon curve can be calculated in nuclear matter if the 
nucleon effective mass is given in the relativistic Hartree approach(RHA) with the inclusion 
of the vacuum fluctuation effects. A brief summary about RHA with density-dependent 
effective coupling constants is given in Sections II and III. A resonable family of the nucleon 
effective mass curves is used to collect the density-dependent coupling constant curves that 
reproduce saturation properties of nuclear matter. These curves define the upper and lower 
bounds of the density-dependent coupling constants in RHA where one-baryon-loop vacuum 
effects are taken into account. The resulting range of the density-dependent coupling con- 
stant would restrict a possible equation of state of a neutron matter or a nuclear system at 
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high density. It also would limit a possible nuclear incompressibility at low density such as 
at the surface of finite nuclei which is important in study of giant resonances related with 
surface vibration. 



II. THE APPLICATION OF THE THERMODYNAMIC LAW IN NUCLEAR 

MATTER 

The nucleon self-energy in nuclear matter has the general form by assuming parity con- 
servation, time- reversal invariance, and hermiticity 0: 

£(Jfe) = £ S (A;, k°, k F ) - 7°S°(i, k°, k F )+j- kX v (k, k°, k F ), (1) 

where k = (k°, k) is the nucleon four momentum and 7 M is Dirac matrix. kp is Fermi 
momentum of nucleon in nuclear matter and k— |k|, so that < k < k F . The self-energy 
can be reduced within the Dirac equation to a potential 

U(k, k F ) = A{k, k F ) + -y°B(k, k F ) (2) 

where the scalar potential A and the vector potential B are given by 

A(k, k F ) = i + ^ v (3) 

jpyv _ yO 

B(k,k F )= i + . (4) 

Here, M is the free nucleon mass and E is the single particle energy of a nucleon in nuclear 
matter including the rest mass energy of a free nucleon. In a mean field approximation, A 
is independent of k and define the nucleon effective mass M*(k F ) as 

M*{k F ) = M + A[k F ). (5) 

The nucleon in nuclear matter behaves like a free particle with energy E* and mass M*, 
where the relation 

E + S° 



E* = E-B= „ ' " = V k 2 + M* 2 (6) 
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holds in the Dirac equation. 

In the thermodynamic treatment, we shall regard nuclear matter as the system contain- 
ing the same kind of particles, namely single-component system. Then, the fundamental 
thermodynamic law is given by 

dE int = TdS - PdV + fidN, (7) 

which describes the change in the internal energy Ei nt arising from small change in the 
entropy S, the volume V, and the number of nucleons N in V. We confine ourselves to 
the ground state nuclear matter here, so the temperature T is zero. Furthermore, the 
Hugenholtz-Van Hove (HV) theorem states that the chemical potential \x equals the highest 
occupied single particle energy Q. Therefore, we have 

— + E A = e(M\k F ) J (8) 
Pb 

where e(M*,k) = E — M is the single particle energy, p B = dN/dV the baryon density and 
E A = dE int /dN the energy per nucleon. 
Using the relation 

p = <> B ^ 0) 

dp B 

and 

Pb = ^k F 3 (10) 

where 7 is the spin-isospin degeneracy, it immediately follows from Eq. (H) that 

E A + ^^ = e(M*,k F ). (11) 
6 dkp 

On the other hand, from Eqs.(f|) and (§) the single particle energy e(M*,k) can be 
expressed as 



e(M*,k) = E — M = \Jk 2 + M* 2 + B — M. (12) 



Substituting Eq. fll2|) for k = kp into Eq. (|TT|) , we have 

k F dEA 



E a + -T-JT 1 = W + M* + B — M. (13) 

O UKp 
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III. THE NUCLEON BINDING ENERGY IN NUCLEAR MATTER 



The role of quantum corrections is crucial in the application of hadronic fields theories 
to nuclear matter. We will use the relativistic Hartree approximation(RHA) which includes 
one-baryon-loop vacuum effects. The RHA involves divergent integrals over the occupied 
negative energy states. The divergences can be rendered finite by including the appropriate 
counterterms in the MFT Largrangian and by defining a set of renormalization conditions. 
The vacuum fluctuation corrections to the binding energy per nucleon can be obtained by 
carrying out the renormalization procedure. The result of the binding energy per nucleon 
with the vacuum corrections is given by @ 



B x , m 2 (M-M*) 2 7 
2 2gj p B {2ir) s pB 



f F d 3 k^Jk 2 + M* 2 + AE VF , (14) 

J 



where 



AE 



VF 



8k F 



M* \ 7 
W* 4 In f —J + M 3 (M - M*) - -M 2 (M - M*) 2 

13 25 

H M(M - M*f (M - M*) 4 

3 X 2 



(15) 



No counterterm for vector self-energy to be finite is needed and therefore the RHA result 
for the vector self-energy term is identical to that of the MFT which is given by 
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For given coupling constants, (g s /m s ) 2 and (g v /m v ) 2 as functions of kp, Eq.(|T^) is a function 
of the effective mass M* and the Fermi momentum kp. The nuclear matter saturation is 
determined by minimizing Ea in kp-M* space. 



For each value of kp, Ea given by Eq. (Jj) should take a minimum at the right value of 
M*, so that the partial derivative of Ea with respect to M* vanishes: 



M -M* 
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m 2 (27r) 3 Jo 



d 3 k 



M* 



k 2 + M* 



m 2 it 2 4 



M* d \n 



M* 



M 2 (M* — M) — -M(M* - M) 2 - — (M* - MY 
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(17) 



By using Eq. ([17]) the factor m 2 / g 2 can be eliminated in Eq. (|T4] ) to have the expression of 
Ea as 

^ B 3 
£ A = — — M + 
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4 + kpMM* - hz F M*' 2 j v'7, 7 ■ M 



l-M 



MM* In 



fcF + y/g + M* 2 
M* 



: M - MM In 



M* 
M~ 



1 3 1 19 1 

•-M 3 (M - M*) + -M 2 (M - M*) 2 + -M(M - M*) 3 - — (M - M*) 4 j 



(18) 



We combine Eq. (|T8| ) with Eq. ([U|) to eliminate £> (i.e., g 2 /m 2 ) and obtain the following 
relation: 



d /£ 4 \ 3M + 3a/H + M* 2 9 



dkF \ kp 



kp 



2,kp 



1 



kp + kpMM* - ^k F M* 2 ) ^ + M 



* 2 



fl -M* 4 - MM* 3 ) In 



fcF + yfc 2 , + m* 2 

M* 



+ 



-M* 4 - MM* 3 ) In 



M* 
M~ 



1 3 1 19 

-M 3 (M - M*) + -M 2 (M - M*) 2 + -M(M - M*) 3 - — (M - M*) 4 



(19) 

One can solve Eq. (|19|) numerically for M*(/c ) for a given value of £^(^0) — —15.75 MeV 
by using the fact that dEA/dkp = at kp = k where k is the saturation density given by 
1.42 fin" 1 . 



IV. CALCULATION OF DENSITY-DEPENDENT COUPLING CONSTANTS 



The scalar meson coupling constant and the nucleon effective mass are related to each 
other through Eq. (|TTD with 7 = 4. The density-dependent coupling constant curve for the 
scalar meson can be obtained immediately when M*(kp) is given: 



g 2 _ tt 2 (M-M*) 
m 2 F(k F ,M*) 



where 



F(kp, M*) 
M* 



M*k F Jkl + M* 2 - M* 3 In 



-M* 3 ln 



kp + ^k 2 F + M* 
M* 
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} + M 2 (M* - M) + -M(M* - M) 2 + — (M* - M) 3 . 



(20) 



(21) 



Eq. (pip shows there is a discontinuity in g 2 s changing from positive infinity to a negative 
infinity at some kp larger than k . This discontinuity results in the change of g s from real to 
a pure imaginary and from attractive n-n scalar coupling to a repulsive n-n scalar coupling. 
To describe a long range nuclear attraction and a short range nuclear repulsion within the 
Walecka's a-u model, the coupling constants g s and g v must be real. If we assume that the 
same property of nuclear interaction, i.e., the attractive scalar interaction and the repulsive 
vector interaction holds at high density too, then g s and g v must stay real. As far as g s is 
real F(kp,M*) should be positive as can be seen in Eq. (|20D . One can limit M*(kp) such 
that F(k F ,M*) > by obtaining the solution of F(k F ,M*) = for M*(k F ) which plays 
a role for the boundary line. The boundary line is shown as the dash-double-dot line in 
Fig. 1. The boundary line is denoted by M'(k F ) as a function of kp. The region below the 
boundary line M'(kp) in Fig. 1 is forbidden for the nucleon effective mass. Eq. ([19]) can be 
integrated numerically to obtain EA{k F ) if M*(k F ) is known. 

We choose the form of M* as a function of kp in the region above the boundary line 
M'(kp) in Fig. 1: 



M*(k F ) 



M 



1 + 



Af-M*(fc ) , k F . n 1 
M*(fe ) 1 fc ; 

1+/3 



kp < k 



[M*(ko) - M'(kp)] + M'(kp), k F > k 



(22) 



where a, f3 and 5 are related by 

1 



+ 



M*(k )-M'(k ) 
{2M*(k ) -M}a 



2M*{k ){M - M*{k Q )}a dM' 

h 2/^0 



M 



dki 



kp=kg 



M 



M*(k ){M - M*(k )}a 



(23) 
(24) 



M{M*(k ) - M'(k )}6 - M*(k ){M - M*(k )}a 

where we choose a > 1 such that dM*/dk F = at k F = 0. We note here that dM*/dk F = 
— oo if a < 1. As discussed above we have chosen the M* curve in the region above the 
boundary line M'(k F ) in the M* — kp plane to confine our calculations to real g s . Eq. ( |22|) 
for M*{k F ) leads to the positive effective scalar coupling constants g1/m 2 s for all a > 1.0 as 



is required. Since g v must also be real, the region of the parameter a (a ^ 3.5) that leads 
to negative g 2 /m 2 , as will be shown, should be discarded. Therefore, for the possible M* 
curve we have chosen Eq. ( p2|) with 1.0 < a < 3.5. The other parameters (3 and 5 have 
been related to a in such a way that M* curve is differentiable at fop = &o up to the second 
order. Eq. (22) has been chosen such that M*(0) = M and M*(/c ) leads to the empirical 



saturation point ( k = 1.42 fm" 1 and E A = -15.75 MeV ). The M* curves are shown for 
each value of a in Fig. 1. 

In numerical calculations, M* given by Eq. (|22|) for a fixed a is substituted in Eq. 
(|i~9l) and the integration with respect to kp is performed to obtain the energy per nucleon 
such that the curve passes through the empirical saturation point. The nucleon binding 
energy curves obtained with M*(/cp) for different a's are plotted in Fig. 2. The calculated 
incompressibility of nuclear matter increases with increasing a (K" 1 = 61, 313 and 565 MeV 
for a = 1.0, 2.25 and 3.5, respectively). 

Once M*(kp) is given by Eq. ( |22"D for a parameter a, the scalar effective coupling constant 
9s/ m 's can be obtained from Eq. (|2"0"|). The g 2 s /m 2 s curves for different a's are shown in Fig. 

3. The effective coupling constant for the vector meson can be expressed in terms of kp, E A 
and M* by using Eqs. (|I(J) and (|T%|). The calculated curves for g^/ml are depicted in Fig. 

4. The negative value of g% at low densities for a > 3.5 gives the upper bound of a. 
Using the density-dependent coupling constants we obtained the equations of state of 

neutron matter as shown in Fig. 5, where these are compared with the equation of state 
of neutron matter with constant coupling constants 0. The saturation density and incom- 
pressibility of neutron matter increase with the increasing a. The calculated results show 
that neutron matter is bound in the regions of 1 < a < 1.65 and 2.79 < a < 3.5. For small 
a, the scalar coupling constant gl/m 2 s and the vector coupling constant g 2 /m 2 are much 
stronger at low density and weaker at high kp than the fixed coupling constants. Thus the 
neutron matter is bound weakly with small kp for small a. On the other hand, for large 
a, both the scalar and vector couplings are weaker than the fixed coupling constants at 
low density (small /cp) and the vector coupling constant g 2 /m 2 is the similar order as the 
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fixed coupling constant at high density while the scalar coupling constant g 2 s /m 2 s becomes 
much stronger than the fixed coupling constant at high density. This density dependence 
of the coupling constants for large a leads the neutron matter to be strongly bound (11.1 
MeV binding energy per neutron). In the mid a value region, the neutron matter is un- 
bound. However, in the region near a = 2.3, a neutron matter has a resonance state at 
&f = ko — 1.42 fm _1 because the coupling constants (g 2 /m 2 , gl/m 2 ) for all the cases (a's) 
are required to be the same at the density. 

V. CONCLUSION 

The results (Figs. 1-4) show that, at low density, the density-dependent coupling 
constants are very sensitive to the choice of the effective mass. However, the nuclear equation 
of state Ea^f) is quite insensitive to the coupling constants at low density. This behavior 
of insensitivity reflects the fact that the inter-nucleon separation is large at low density. 
Due to the requirement of the empirical saturation in our calculation, all the choices of the 
effective mass, of course, give the same energy and coupling constants at the saturation 
point. However the nuclear compressibilities are all different. 

We can see from the results that the equation of state is quite sensitive to the coupling 
constants at high density as we expected. In addition, the results show that the equation 
of state of a nuclear system is also sensitive to the values of the coupling constants at the 
intermediate density, around half the saturation density. This behavior suggests us that the 
detailed study of surface vibration modes of finite nuclei is very important in extrapolating 
the nuclear equation of state to a high density. To pursue a research in this direction, we 
need futher investigation of our starting point of calculation whereas here we started with 
the effective mass functions M*(kp) just for simplicity. The minimum requirement in this 
search is that all the fits should pass through at kp = and kp = ko, i.e., a zero density 
limit and a nuclear saturation point. 

Fig. 5 also shows that the density-dependent coupling constants allow the existence of 
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neutron matter within RHA. The equations of state of neutron matter with the density- 
dependent coupling constants show saturation points with different incompressibilities of 
neutron matter. The saturation points for different a's are compared with that for constant 
coupling constants || in Fig. 5. The saturation state appears at higher density and the in- 
compressibility is higher if a takes on the larger value. We note here that both the resonance 
states with positive saturation energy and the bound states with negative saturation energy 
appear in the neutron matter calculations with the density-dependent coupling constants. 
The saturation energies are positive in the coupling constant region 1.65 < a < 2.79 where 
RHA does not allow a stable neutron star. However, for the regions 1.0 < a < 1.65 and 
2.79 < a < 3.5 RHA provides bound states of neutron matter, and for a = 3.5 the binding 
energy appears to be 11.1 MeV the incompressibility to be 1258 MeV at hp = 1.94 fm _1 . 
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Research Institute Grant No. 95-2422. One of the authors (S.J.L.) would also like to 
acknowledge the Korea Science and Engineering Foundation for financial support under 
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FIGURES 

FIG. 1. Nucleon effective masses M* as functions of Fermi momenta &f for nuclear matter. 



The M* curves are chosen as Eq. (22) with a parameter a such that they lead to the empirical 
saturation point in nuclear matter. 

FIG. 2. Binding energies per nucleon in nuclear matter calculated from the nucleon effective 
masses given by Eq. (22). 

FIG. 3. Density-dependent effective coupling constant curves for the cr-meson exchange. The 
constants g1/m 2 s are obtained by using Eq. (|20|) with the nucleon effective masses given by Eq. 



FIG. 4. Density-dependent effective coupling constant curves for the u;-meson exchange. 

FIG. 5. Binding energies per nucleon for neutron matter (7 = 2) calculated from the den- 
sity-dependent coupling constants gl/m 2 in Fig. 3 and g 2 /m 2 in Fig. 4. The maximum binding 
energy per neutron is 11.1 MeV for a = 3.5. 



12 



